Inelastic neutron scattering study of magnetic excitations in Sr 2 Ru04. 
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Magnetic excitations in Sr2Ru04 have been studied by inelastic neutron scattering. The magnetic 
fluctuations are dominated by incommensurate peaks related to the Fermi surface nesting of the 
quasi-one-dimensional d xz - and d yz -bands. The shape of the incommensurate signal agrees well with 
RPA calculations. At the incommensurate Q-positions the energy spectrum considerably softens 
upon cooling pointing to a close magnetic instability : Sr2Ru04 does not exhibit quantum criticality 
but is very close to it. w/T-scaling may be fitted to the data for temperatures above 30 K. Below 
the superconducting transition, the magnetic response at the nesting signal is not found to change 
in the energy range down to 0.4meV. 

PACS numbers: PACS numbers: 78.70.Nx, 75.40.Gb, 74.70.-b 



I. INTRODUCTION 

Sr 2 Ru04 is still the only superconducting layered per- 
ovskite besides the cupratesQ; however, in contrast to 
the cuprate high temperature superconductors (HTSC), 
superconductivity in Sr2Ru04 develops in a well defined 
Fermi-liquid state 0, ||] . Nevertheless the superconduct- 
ing state and the pairing mechanism in S^RuC^ are 
highly unconventional. The present interest in this com- 
pound goes far beyond the simple comparison with the 
cuprate high temperature superconductors. 

The extreme sensitivity of the superconducting tran- 
sition temperature on non-magnetic impurities suggests 
a non-phonon mechanism It is further established 
that superconductivity in Si~2Ru04 is made of spin- 
triplet Cooper pairs and breaks time-reversal symmetry 

H . The strongest experimental argument for that can 
be found in the uniform susceptibility measured either 
by the NMR-Knight-shift or polarized neutrons experi- 
ments j| | and in the appearance of spontaneous fields 
in the superconducting state reported by fiSK JtJ. A 
spin-triplet p-wave order parameter had been proposed 
before these experiments ||, in the idea that supercon- 
ductivity arrises from a dominant interaction with ferro- 
magnetic fluctuations in analogy to superfluid 3 He. Rice 
and Sigrist stressed the comparison with the perovskites 
SrRu03 and CaRu03 which order ferromagnetically or 
are nearly ferromagnetic respectively 0. Evidence for 
ferromagnetic fluctuations in S^RuC^ was inferred from 
NMR-experiments : Imai et al. found that Q f the O- 
and of the Ru-NMR exhibit a similar temperature de- 
pendence and interpreted that this could be only due to 
ferromagnetic fluctuations jlO| . 

The macroscopic susceptibility in Sr2Ru04 is en- 
hanced when compared with the band structure calcu- 
lation but only weakly, in particular its temperature 
dependence remains flat §, [Uj]. There exist also lay- 



ered ruthenates which are very close to ferromagnetic 
order at low temperatures : the purest Sr3Ru207 sam- 
ples show meta-magnetism and samples with some- 
how less quality even order ferromagnetically p2[ . A 
highly enhanced susceptibility pointing to a ferromag- 
netic instability is also observed in the phase diagram of 
Ca2- a; Sr a ;Ru04 Q but for a rather high Ca con- 
centration, Cai.5Sro.5Ru04 (IB). In these nearly ordered 
layered ruthenates, the susceptibility is about two orders 
of magnitude higher than that in Sr2Ru04 and strongly 
temperature dependent. 

Some doubt about the unique role of ferromagnetism 
in Sr2Ru04 arose from the strong moment antiferro- 
magnetic order observed in the Ca analog ]l6| ], which 
inspired Mazin and Singh to perform a calculation of the 
generalized susceptibility based on the electronic band 
structure fll7| . Surprisingly they found that the domi- 
nating part is neither ferro- nor antifcrromagnetic but 
incommensurate. The Fermi-surface in Sr2Ru04 is well 
studied both by experiment |L8| and by theory jl9| with 
satisfactory agreement. Three bands are contributing 
to the Fermi-surface which may be roughly attributed 
to the three t2 S -levels, the d xy -, d yz - and d^-orbitals, 
occupied by the four 4d-electrons of Ru 4+ . The d xy - 
orbitals hybridize in the xy-plane and, therefore, form a 
band with two-dimensional character, 7-band. In con- 
trast, the d xz and the d yz orbitals may hybridize only 
along the x and the y-directions, respectively, and form 
quasi-one-dimensional bands, a-band and /3-band, with 
flat sheets in the Fermi-surface, a-sheet and /3-sheet. 
The latter give rise to strong nesting and enhanced 
susceptibility for q = (0.33,^,0) or q = (q-r, 0.33,0) 
|]l7|| . Along the diagonal both effects strengthen each 
other yielding a pronounced peak in the susceptibility at 
qi-cai = (0.33,0.33,0). Using inelastic neutron scatter- 
ing (INS) we have perfectly confirmed this nesting sce- 
nario p0| . The dynamic susceptibility at moderate ener- 
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gies is indeed dominated by the incommensurate fluctu- 
ations occurring very close to the calculated position at 
qi = (0.30, 0.30, 0). Furthermore, we found a pronounced 
temperature dependence which can explain the temper- 
ature dependence of both Ru- and O- jry-NMR |Tc| ]. 

In the meanwhile, several groups have worked on theo- 
ries where superconductivity is related to the incommen- 
surate fluctuations ^2], ^3|. These theories, however, 
do not yet give an explanation for the fine structure of the 
order parameter. Specific heat data on the highest qual- 
ity samples clearly indicates the presence of line nodes in 
the superconducting state (24) . Ultrasonic and ther- 
mal conductivity ]26| results would disagree with vertical 
line nodes which leave horizontal line nodes as the only 
possible ones. Zhitomirsky and Rice [^7j assume that su- 
perconductivity is transferred from the 7-band to the 
a-band and /3-band by a proximity effect and get a con- 
clusive explanation for the horizontal line nodes. In this 
model superconductivity should be mainly related to ex- 
citations associated with the active 7-band, which so far 
have not been characterized. Therefore, it appears still 
interesting to further deepen the study of the magnetic 
fluctuations in Sr2Ru04 . 

In this paper we report on additional INS studies in 
Sr2RuC>4 in the normal as well as in the superconducting 
phase. We present a more quantitative analysis of the 
incommensurate fluctuations related to the a-sheet and 
/3-sheet and discuss the possible contributions due to the 
two-dimensional 7-band. 



suppress higher order contamination, the final neutron 
energy was fixed at 14.7meV. All diaphragms determin- 
ing the beam paths were opened more widely than usu- 
ally in order to relax Q-resolution, since the magnetic 
signals are not sharp in Q-space. In most experiments, 
the scattering plane was defined by (1,0,0) and (0,1,0) 
directions in order to span any directions within the (ab) 
plane. An additional experiment has been made with the 
(1,1,0) and (0,0,1) directions with the plane to follow the 
spin fluctuations along the c* axis. 

Studies aiming at the response in the magnetic exci- 
tation spectrum on the opening of the superconducting 
gap require an energy resolution better than the expected 
value for twice the superconducting gap. Therefore, such 
experiment is better placed on a cold triple-axis spec- 
trometer even though this implies a sensitive reduction 
in the flux. We have made preliminary studies using 
the cold spectrometers 4F at the Orphee reactor and a 
more extensive investigation on the spectrometer IN14 
at the ILL. These instruments possess PG monochroma- 
tors (double at 4F and focusing at IN14) and focusing 
analyzers. The final neutron energy was fixed at Ef = 5 
meV on both cold source spectrometers where a Be filter 
has been employed to cut down higher-order wavelength 
neutrons. Cooling was achieved by use of a dilution- and 
a He 3 -insert at 4F and IN14 respectively. 



B. Theoretical background for magnetic neutron 
scattering 



II. EXPERIMENTAL 

A. Experimental setup 

Single crystals of Sr2Ru04 were grown by a float- 
ing zone method in an image furnace; they exhibit the 
superconducting transition at T c =0.7, 1.35 and 1.43K 
and have volumes of about 450mm 3 each. Since most of 
the measurements were performed in the normal phase, 
where the differences in T c should not affect the mag- 
netic excitation spectrum, we aligned the three crystals 
together in order to gain counting statistics in the INS 
experiments. Count rates in the ruthenates experiments 
are relatively small already due to the steeper decrease 
of the Ru magnetic form- factor with increasing scattering 
vector. For the measurements below and across the su- 
perconducting transition we mounted only the two crys- 
tals with relatively high T c together. The mounting of 
the two or three crystals was achieved with individual 
goniometers yielding a total mosaic spread below 0.5 de- 
grees. 

We used the thermal triple axis spectrometer IT in- 
stalled at the Orphee-reactor (Saclay, France) in order to 
further characterize the scattering in the normal state. 
The instrument was operated with double focusing py- 
rolithic graphite (PG) monochromator and analyzer, in 
addition PG filters in front of the analyzer were used to 



The magnetic neutron cross section per formula unit 
is written in terms of the Fourier transform of the spin 
correlation function, S a/ 3(Q, oj) (labels a, (3 correspond to 
x,y,z) as H; 
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where k^ and k/ are the incident and final neutron wave 
vectors, r§=0.292 barn, F(Q) is the magnetic form fac- 
tor. The scattering-vector Q can be split into Q =q+G, 
where q lies in the first Brillouin-zone and G is a zone- 
center. All reciprocal space coordinates (Q x , Qy, Qz) a- re 
given in reduced lattice units of 27r/a or 2-k/c. 

According to the fluctuation-dissipation theorem, the 
spin correlation function is related to the imaginary part 
of the dynamical magnetic susceptibility times the by one 
enhanced Bose-factor: 
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In case of weak anisotropy, which is usually observed in a 
paramagnetic state, Xap" (Q> reduces to x" (Q> ^)5 a ,i3- 
Note that for itinerant magnets, anisotropy of the suscep- 
tibility can occur due to spin-orbit coupling. F(Q) can 
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be described by the Ru + magnetic form factor in first ap- 
proximation. Once determined, the magnetic response is 
converted to the dynamical susceptibility and calibrated 
in absolute units by comparison with acoustic phonons, 
using a standard procedure depicted in p0[. 



III. RESULTS AND DISCUSSION 

A. RPA analysis of the magnetic excitations 

At low temperature Sr2Ru04 exhibits well defined 
Fermi-liquid properties; therefore, it seems appropriate 
to assess its magnetic excitations within a RPA approach 
basing on the band structure. Density functional calcula- 
tions in LDA were performed by several groups and yield 
good agreement with the Fermi-surface determined in de- 
Haas van Alphen experiments. The bare non- interacting 
susceptibility, x°(q), can be obtained by summing the 
matrix elements for an electron hole excitation [E8fl: 



X°(q,w) = -J2 M ^- 



/(£k+ qj ) - /Ok,i) 



£k-| 



IC 



(3) 



where e — >0, / is the Fermi distribution function and 
Ek the quasiparticle dispersion relation. This was first 
calculated by Mazin and Singh [fL7| under the assumption 
that only excitations within the same orbital character 
are relevant (the matrix-elements Mi^k+qj are equal to 
one or zero) . Mazin and Singh predicted the existence of 
peaks in the real part of the bare susceptibility at lu=0 
due to the pronounced nesting of the a-bands and fi- 
bands. These peaks were calculated at (0.33,0.33,0) and 
experimentally confirmed very close to this position at 
qi=(0.3,0.3,0), see Fig. 1. In addition to the peaks at qi, 
this study finds ridges of high susceptibility at (0.3, q y ,0) 
for 0.3 < q y < 0.5 and some shoulder for < q y < 0.3. 

The susceptibility gets enhanced through the Stoner- 
like interaction which is treated in RPA by : 



X(q) 



X°(q) 



l-7(q)x°(q) 



(4) 



with the q-dependent interaction /(q). For the nesting 
positions Mazin and Singh get /(q)x°(q)=l-02, which 
already implies a diverging susceptibility and a magnetic 
instability. 

In Fig. 1 we show a scheme of the (hkO)-plane in re- 
ciprocal space. Due to the body centering in space-group 
I4/mmm any (hkl)-Bragg-points have to fulfill the condi- 
tion (h+k+1) even; (100) for instance is not a zone-center 
but a Z-point. The boundaries of the Brillouin-zones are 
included in the figure. Large filled circles indicate the 
position of the incommensurate peaks as predicted by 
the band-structure and as observed in our previous work. 
The thick lines connecting four of them correspond to the 
ridges of enhanced susceptibility also suggested in refer- 
ence |Et|. 
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FIG. 1: Upper part : Scheme of the (hkO)-plane in reciprocal 
space. Thin lines show the boundaries of the Brillouin-zones 
and small filled and open circles the zone-centers and the Z- 
points of the body centered lattice. Large filled circles indi- 
cate the position of the incommensurate peaks and thick lines 
connecting four of them correspond to the walls of enhanced 
susceptibility also suggested in reference [17]. The dashed 
double arrows show the paths of the constant energy scans 
frequently used in this work: along [010]-direction - y-scan, 
transverse to the Q- vector, r-scan, and along the diagonals in 
[110]-direction, d-scans. Lower part : Imaginary part of the 
generalized susceptibility calculated by RPA along the three 
paths indicated in the upper part as dotted lines. 



In the meanwhile several groups have performed simi- 
lar calculations which all agree on the dominant incom- 
mensurate fluctuations |23, 29, 30|, |3l|, |32[|. However there 
are serious discrepancies concerning the detailed struc- 
ture of the spin susceptibility away from q;. These dis- 
crepancies mainly rely on the parameters used to describe 
the electronic band structure, on the choice of I(q) and 
on the inclusion of more subtle effects such as spin-orbit 
or Hund couplings. 

In order to compare directly to the INS experiments, 
see equations (1) and (2), it is necessary to perform the 
RPA analysis by taking into account both real and imag- 
inary part of the susceptibility. Morr et al. |}2| report 
such calculations obtained by fitting the band structure 
to the ARPES results p3) . They find in addition to the 
peak at qi quite strong intensity near Pi=(0.3,0.5,0), 
(the middle of the walls, see Fig. 1), which is even com- 
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parable to that at qi in the bare susceptibility. Simi- 
lar results were obtained in references |2^, |2{|. Ng and 
Sigrist [|o| find much less spectral weight in the ridges 
at (0.3, q y ) for 0.3 < q y < 0.5 but stronger shoulders 
< q y < 0.3. In addition, they calculate the separate 
contribution of the 7-band which does not show a par- 
ticular enhancement in the ferromagnetic q-range but is 
little structured. Eremin et al. |3l|| calculate the sus- 
ceptibility taking into account strong hybridization and 
obtain results somehow differnt from the other groups. 
They find a strong signal at in addition there is some 
enhancement of the susceptibility related to the van-Hove 
singularity of the 7-band. This contribution occurs quite 
close to the zone center at q V H=(0.15,0,0). 

We have performed the full RPA analysis basing on 
the LDA band structure reported in reference jl7| [n]] 
in order to accompany our experimental investigations. 
We first calculate the bare electron hole susceptibility 
X° from the usual expression, equation (3). For the band 
energies e(k x ,k y ,k z ) we use the expressions of Mazin and 
Singh Jl9[ for the three mutually non hybridizing tight- 
binding bands in the vicinity of the Fermi level: 

e xy (k x , k y , k z ) = A00mev(—1 + 2(cos(k x ) + cos(k y )) 
— 1.2cos{k x )cos(ky)) (5) 

e xz (k x , k y , k z ) — 400mev(— .75 + 1.25cos(fe x ) — 

.5cos(k x /2)cos(k y /2)cos(k z /2)) (6) 

£yz{k x ,ky,k z ) = 400meu(— .75 + 1.25cos(k y ) — 

.5cos(k x /2)cos(k y /2)cos(k z /2)) (7) 

We also used the crude approximation that matrix- 
elements for transitions between bands of the same char- 
acter are equal to one and others zero. The q-dependent 
"ferromagnetic" interaction (Stoner factor) I(q) is taken 
to be equal to (following reference 0): I(q) = fqi^jj^j?- 
With this choice of I(q) the calculated static susceptibil- 
ity x(q = 0, ui = 0) is slightly lower than the measured 
one. If I(q) is chosen larger an instability appears at the 
incommensurate wave vector. 

Our results for the imaginary part of the generalized 
susceptibility at an energy transfer of 6meV are given 
in the lower part of figure 1. Besides the dominating 
nesting peak near qi there is a further contribution near 
(0.15,0.15,0) which is related to the 7-sheet. In contrast 
we find a small susceptibility near Pi and for (q x fl,0) 
with small values of q x . 

Since the q-position of the magnetic excitations were 
found not to depend on energy, it is easiest to observe the 
signal in INS by scanning at constant energy. The scan 
paths are included in Fig. 1, they are purely transverse 
or rocking-like, r-scan, along a [100]-direction, x-scans, 
and in diagonal direction, d-scans. 

The observed signal is rather broad and, therefore, the 
scans performed are extremely wide covering complete 
cuts through the Brillouin-zone. This further implies 
that the background (BG) may be non-constant at least 
sloping. Also, the signals are relatively weak compared 
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to typical triple axis spectrometry problems; this implies 
that sample independent BG-contributions which usually 
are negligible play a role. Fig. 2 presents the results of 
d-scans at different temperatures clearly demonstrating 
the gain in statistics compared to the previous work p0[ . 
The magnetic intensity shown in Fig. 2 disappears upon 
heating but this effect gets partially compensated by the 
gain through the Bose-factor. 



B. Shape of the incommensurate signal 

The fact that the incommensurate signal around a 
zone-center and around a Z-point, (001), are equivalent 
already indicates that the coupling between Ru02-plancs 
is negligible, i.e. that in-phase and out-of-phase cou- 
pling are indistinguishable. The 2D-character has been 
directly documented by Servant et al. J34| who found 
no q;-dependence at (0.3,0.3,(7;) f° r qz between -0.5 and 
0.5. We find the same result by varying Q; in a broader 
range between 2 and 5 in (0.3,0.3,Q;), see Fig. 3. This 
2D character is actually surprising since the dispersion 
relation of the quasi-particles involved in the computa- 
tion of x"(q, w) is not purely 2D fl7fl . One may therefore 
expect weak spin correlation along the c*-direction. 

Recently, it has been shown that these fluctuations 
freeze out into a spin density wave (SDW) ordering by 
a minor replacement of Ru by Ti j35j. In this order- 
ing a very short correlation length along the c-direction 
nicely reflects the 2D-character of the incommensurate 
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FIG. 3: Height of the magnetic signal at the incommensurate 
position as function of Qi; the dashed line shows the depen- 
dence expected from the Ru + -form factor and the solid line 
that assuming the spin-density distribution observed recently 
in Cai.sSro.sRuCUjSSl , see text. 



inelastic signal. The SDW propagation vector finally fa- 
vored in the static ordering corresponds to the out-of- 
phase coupling between neighboring layers. This static 
interlayer-coupling might be even due to the CDW always 
coupled to a SDW, which, however, has not yet been dis- 
covered so far. One may add that the stripe ordering in 
Lai.67Sro.33Ni04, which differs from the Si^Rui-zTi^C^- 
case in many aspects but nevertheless may be consid- 
ered as a mixed SDW-CDW-ordering, occurs at the same 
propagation vector J3(| . 

The wider Qj-dependence of the incommensurate sig- 
nal shown in Fig. 3 can give information about the 
anisotropy of the excitation, since INS measures only 
the spin component perpendicular to Q. We consider the 
diagonal susceptibility x" a /3 with tetragonal symmetry 
X± ■= X"xx = x"yy X" zz- The measured intensity is 
then given by : 
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which implies that the observations at high Q/ favor the 
in-plane component of the susceptibility. For the de- 
tailed analysis, one has to compare with the form fac- 
tor. In the figure we show the Qj-dependence assuming 
that spin density is localized at the Ru-site and may be 
modeled by the form factor of Ru 1+ (3?]] . The Ru-form- 
factor dependence underestimates the signal at higher 
Qz-values; however, the Ru 1+ -form factor is certainly a 
too crude approximation. Measurements of the spin- 
density distribution induced by an external field have not 
been very precise due to the small magnetic susceptibility 
and the resulting small moment in Sr2Ru04 |3q| . How- 
ever, in Cai.5Sro.5Ru04 which exhibits ferromagnetic 
ordering below IK and whose low temperature suscepti- 
bility is about two orders of magnitude higher than that 
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denote fits with Gaussians and dashed lines the calculated 
imaginary part of the generalized susceptibility at the energy 
of 4.1meV (shifted in y-axis). Arrows indicate the position of 
the shoulder seen in the RPA calculations. 



in Sr2RuC>4 |T^, 14, |T^], it has been possible to study 
the field induced spin-density distribution. These exper- 
iments revealed an extremely high amount of spin-density 
at the oxygen position, about one third of the total mo- 
ment pfj| and an orbital contribution at the Ru-site. By 
use of the Cai.5Sro.5Ru04 spin-density distribution one 
obtains a good description of the Q/-dependence given in 
Fig. 3. However, the form-factor in Sr2Ru04 should 
be even more complex. Since the main contribution 
originates from the flat e^y-orbitals there must be 
an anisotropy in the effective form factor which indeed 
was observed in Cai.sSrp.sRuC^flfjt . Qualitatively, the 
anisotropy in the form-factor has to be compensated by 
some weak anisotropy in the spin susceptibility, i.e. by 
an enhanced out-of-plane component. Such susceptibility 
anisotropy corresponds to the orientation of the spins in 
the SDW ordering phase in Sr2Rui_ x Ti x 04 Q where 
the spins are aligned parallel to the c-direction and also 
to the conclusion deduced by Ng and Sigrist from the 
influence of spin-orbit coupling in Sr2RuO4[|30|. 

- Shoulder of the incommensurate peak - In order to 
examine the shape of the incommensurate peak and to 
verify the existence of the ridges of additional nesting in- 
tensity or the additional peak Pi, reported in different 
band structure analyzes |29[ [nj |3^| , we have 
scanned across qi along the [100] or [010]-directions, 
y-scans see Fig. 1. The results are shown in Fig. 4. 
One can recognize that the incommensurate peak is not 
symmetric but exhibits always a shoulder to the lower 
(fc-side in absolute units. The shoulder is seen in many 
scans in reversed focusing conditions of the spectrometer 
configuration excluding an experimental artefact. Our 
full RPA calculation nicely agrees with such shoulder; 
the thin line in Fig. 4 shows the calculated imaginary 
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FIG. 5: Mapping of the INS intensity by constant energy 
scans with fixed Q y at 4.1meV energy transfer. The colour 
plot was obtained by adding symmetrical data after subtrac- 
tion of the scattering angle dependent background. 



part of the susceptibility at the energy transfer of the 
scan and describes the observed signal perfectly besides 
a minor offset in the position of the incommensurate sig- 
nal. In contrast, neither the experiment nor our full RPA 
analysis yield significant intensity in the ridges, i.e. the 
range (g^, 0.3,0) with 0.3 < q x < 0.5, see Fig. 1 and 4. 
In particular an intensity at Pi only three times weaker 
than that at qi would have been easily detected exper- 
imentally. The nesting peak appears to be isolated with 
two shoulders along [100] and [010] to the lower absolute 
cfc^-sides. These shoulders are connected to the 7-sheet. 



C. Possible additional magnetic excitations 

In none of the band-structure calculations there is evi- 
dence for a strong and sharp enhancement of any suscep- 
tibility exactly at the zonc-ccntcrpointing to a ferromag- 
netic instability [jl?], |2^, ||(], [n], |3j. However, several of 
these calculations find some large susceptibility near the 
zone-center which can be associated with the van Hovc- 
singularity of the 7-band closely above the Fermi-level. 
Eremin et al. (3^] report this signal at q=(0.1-0.2,0,0), 
other groups find a small peak along the diagonal (q,q,0) 
p3| , [29| , |30|,. Our own analysis also yields such a signal 
which is found to be strongest at (0.15,0.15,0) and which 
levels out along the [100] and [010] directions, see Fig. 1. 

In order to address this problem we have mapped out 
the intensity for (Q x ,Q y fi) with —0.5 < Q x < —0.0 and 
0.6 < Q y < 1.5; these scans are shown in Fig. 5 af- 
ter subtraction of the scattering angle dependent back- 
ground. It is obvious that the incommensurate peak is 
by far the strongest signal. At low temperature one may 
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intensity (corresponding to the imaginary part of the suscep- 
tibility multiplied with the thermal factor) at different tem- 
peratures and energies, c) sum of the four scans shown in 
part a) for different temperatures at tj=4.1meV. c): Sum of 
the four scans shown in part a) at w=10.3meV (at 295K not 
all the four scans have been performed). . 



roughly estimate any additional signal to be at least a 
factor 6 smaller than the nesting peak. The analysis of 
such weak contributions is quite delicate and demands a 
reliable subtraction of the background. Nevertheless the 
scans in the Fig. 5 indicate some scattering closer to the 
zone center. However, this contribution is not sharply 
peaked at the T-point but forms a broad square or a cir- 
cle. The signal roughly agrees with the prediction that 
the 7-band yields magnetic excitations near the zone- 
center. 

Since the sloping background is a major obstacle to an- 
alyze the additional contributions we tried to compensate 
these effects by scanning from Q=(0 1 0) along the four 
diagonals, which are illustrated in Fig. 1, and by sum- 
ming the four scans. The results are given in Fig. 6. In 
the four single scans, Fig. 6a), one recognizes the nesting 
signal with its intensity being determined by the form- 
factor and the sloping background. The summed scans, 
see Fig. 6c), should have constant background and the 
low temperature summed scan once more documents that 
the nesting signal is by far the strongest one. However, 
upon heating additional scattering contributions seem to 
become enhanced in intensity in particular compared to 
the nesting signal which decreases; note that the Bosc- 
factor will already strengthen any signal by a factor three 
in the data in Fig. 6c). Also, at higher energy the addi- 
tional scattering seems to be stronger as seen in Fig. 6d) 
(the background is strongly sloping even in the sum due 
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to a smaller scattering angle). The energy and temper- 
ature dependence of the additional contribution corre- 
sponds to that predicted by the full RPA analysis for the 
7-band magnetic contribution shown in Fig. 6b). The 
spectral weight at Q=(0. 15,0. 15,0) relative to that of the 
nesting feature increases upon increasing temperature or 
energy as it is expected for a signal directly related to 
the van Hove singularity. Due to the agreement between 
the INS results and the RPA calculations, we suggest to 
interpret the additional broad scattering as being mag- 
netic in origin; however, a polarized neutron study would 
be highly desirable. This scattering further might be rel- 
evant for a quantitative explanation of the NMR-data 
]To| , po| ], in particular its temperature independent part. 

Sr 2 RuC>4 is not close to ferromagnetic order but sub- 
stitution of Sr through Ca yie lds such order for the con- 
centration Cai.5Sr .5RuO4[|15|. This doping effect was 
explained in a band structure calculation ]40[ | as arising 
from a down shift in energy of the 7-band pushing the 
van Hove-singularity closer to the Fermi-level. In such 
samples the 7-band magnetic scattering should there- 
fore become strongly enhanced. Indeed first INS studies 
on these compounds reveal broad signals similar to the 
additional scattering described above, but much stronger 
pl[ . This strongly supports a magnetic interpretation of 
the scattering in figures 5 and 6. 



D. Combined temperature and energy dependence 
of the incommensurate signal 



The energy and temperature dependence has been 
studied in more detail by performing r-scans across the 
incommensurate position, since in this mode the back- 
ground is almost flat. However, due to phonon contribu- 
tions we could not extend the measurements to energies 
higher than 12meV, which is slightly below the lowest 
phonon frequency observed at qi |Q . In particular the 
large scans required to cover the broad magnetic signal 
prevent any analysis within the phonon band frequency 
range on a non-polarized thermal triple axis spectrome- 
ter. 

The results of the scans are given in Fig. 7. At low 
temperature we find an energy spectrum in good agree- 
ment to that published earlier |2C[| . In the range up to 
12mcV we observe at all temperatures an energy inde- 
pendent peak width, which, however, increases upon in- 
crease of temperature. For temperatures much higher 
than 160K the background considerably increases, see 
Fig. 2, and prevents a detailed analysis within reason- 
able beam-time. In Fig. 7b) we show the temperature 
dependence of the peak width averaged over the different 
energies which agrees well to the results obtained from 
the single scans with less statistics reported in reference 
p0| . Even at the lowest temperature the width of the 
signal remains finite. 

The spectral functions have been fitted by a single re- 
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laxor behavior fl43l| 



X"(qi,w) = x'(<u,o)- 



(9) 



where T is the characteristic energy and x'(qi, 0) the am- 
plitude which corresponds to the real part of the gener- 
alized susceptibility at w=0 according to the Kramers- 
Kronig relation. The Q-dependence of the signal may be 
described by a Lorentzian distribution with half width 
at half maximum k, but due to experimental broaden- 
ing the constant energy scans are equally well described 
by Gaussians. At low temperature the spectrum clearly 
exhibits a characteristic energy as seen in the maximum 
of the energy dependence, but this maximum is shifted 
to higher values upon heating. Fig. 7c) and 7d) re- 
port the results of the least squares fits with the single 
relaxor. Although the statistics is still limited, the ten- 
dencies can be obtained unambiguously. The height of 
the spectral functions, x'(QiiO), rapidly decreases upon 
heating, which overcompensates the broadening of the 
signal in Q-space. 
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Our finding that the characteristic energy in the range 
6-9meV is well defined only at low temperatures can be 
related with the far-infrared c-axis reflectance study by 
Hildebrand et al. since the optical spectrum shows 
a resonance in_this energy range at low temperatures. 
In reference 



20 



we have compared the temperature de- 
pendence of the incommensurate signal with that of the 
spin-lattice relaxation rate T\ measured by both 17 and 
101 Ru NMR experiments (hJ. These NMR-techniques 
probe the low energy spin fluctuations (u> — > with re- 
spect to INS measurements); furthermore, they integrate 
the fluctuations in q-space. (1/T\T) is related to the 
generalized susceptibility and the INS results by p5| : 



(1/TiT) 



I^(q)l 



X"(q,w) 



(10) 



with | j4(q) | the hyper fine fields. {1/T\T) corresponds 
hence to the slope of the spectral function in Fig. 7a) 
times the extension of the signal in Q-space. The new 
data perfectly agrees to the former result, the loss of 
the incommensurate signal upon heating may explain al- 
most entirely the temperature dependent contribution to 
(1/TxT) 0. 

The temperature dependence of the magnetic excita- 
tion spectrum at the incommensurate position may be 
analyzed within the results of the self consistent renor- 
malization theory described in reference J43| . In an anti- 
ferromagnetic metal the transition is governed by a sin- 
gle parameter related to the Stoner-enhancement at the 
ordering wave- vector, 8 = 1 — /(qi)x°(qO- The charac- 
teristic entities of the magnetic excitations are then given 
by: 



K 2 OC S 

r oc s 



(11) 

(12) 



in q-space implying a short correlation length of just 3- 
4 lattice spacings. The temperature dependence of the 
magnetic excitations corroborate our recent finding that 
only a small amount of Ti is sufficient to induce SDW 
magnetic ordering 

Since Sr2Ru04 is close to a quantum critical 
point it is tempting to analyze whether the excita- 
tion spectrum is governed by some ^-scaling, as it has 
been claimed for the hi gh temperature superconductor 
La2_ x Sr x Cu0 4 @ @ and for CeCu 5 . 9 Au .i |§. 
One would expect that the susceptibility is given by : 



X"( qi , W ,T)ocT-^(-) 



(14) 



In Fig. 8 we plot the w"x" (qi, w)-data of Fig. 7 and 
that obtained previously as a function of temperature p3] 
against S for a=0.75 and 1.0. Only the data at higher 
temperatures agree with the scaling concept, demonstrat- 
ing that Si'2Ru04 is not a quantum critical point. The 
schematic inset may illustrate the phase diagram, where 
the magnetic transition is determined by some parame- 
ter r (external pressure or composition). At the critical 
transition one would observe quantum criticality in the 
entire temperature range, whereas for r-values where the 
transition is suppressed quantum criticality is observed 
only at higher temperatures. One then may expect a 
cross-over temperature T* where the system transforms 
from an unconventional metal at high temperatures to- 
wards a Fermi-liquid at low temperatures |50(| . Only in 
the temperature range above T* the magnetic excitations 
should exhibit the related ^-scaling. Our data clearly 
shows that such scaling can be fitted to the data only for 
the three higher temperatures studied. The description 
with the scaling concept seems to be slightly better for 
the exponent a=0.75. The temperature dependent data 
suggests a crossover near 30K. This cross-over agrees very 
well to that seen in electronic transport properties, where 
well defined Fermi-liquid behavior is only observed below 
about 25K |,||. 



- 1 r oc 8 . (13) 

When the system approaches the phase transition, the 
unique parameter 8 diminishes, which behavior should 
be observable in all three parameters. Equations (11- 
13) imply a sharpening of the magnetic response in q- 
space as well as in energy and a divergence of the sus- 
ceptibility at the ordering vector. Fig. 7 qualitatively 
confirms this picture. All the relevant parameters, see 
the right scales in Fig. 7, decrease towards low temper- 
ature. Therefore, one may conclude that Si~2Ru04 is 
approaching the SDW transition related to the nesting 
effects upon cooling. However, all these parameters do 
not vanish completely but remain finite even at the low- 
est temperatures in agreement with the well known fact 
that Sr2Ru04 does not exhibit magnetic ordering. In 
particular the magnetic scattering remains rather broad 



E. Magnetic scattering in the superconducting 
phase 

As emphasized by Joynt and Rice |pl| , the wave- 
vector- and energy-dependent spin susceptibility in su- 
perconductors reflects directly the vector structure of the 
superconducting (SC) gap function, allowing a complete 
identification of the SC order parameter symmetry. In- 
elastic neutron scattering experiments have the potential, 
in principle, to determine the superconducting order pa- 
rameter. In high-T c superconductors, the spin-singlet d- 
wave symmetry SC gap induces a striking modification of 
the spin susceptibility in the superconducting state. As 
a consequence, the so-called "magnetic resonance peak" 
has been observed in the superconducting state of various 
copper oxide superconductors by INS[j52], ^3). Therefore, 
a similar experiment in Sr2Ru04WOuld certainly be in- 
structive about the SC gap symmetry. 
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The enhanced spin susceptibility has been calculated 
in H^, [EJ, H|, considering a spin-triplet p-wave super- 
conductiviting state with d(k)=z(k x ± ik y ). Note that 
in such a case, the superconducting gap is isotropic due 
to the particular shape of the Fermi-surface in S^RuOzi. 
For the wavevector qi, Kee et al. |34| and Morr et al. 
p2[ predicted that below T c spectral weight is shifted 
from below twice the superconducting gap, A, into a res- 
onance like feature close to 2A. Morr et al. find the res- 
onance in the zz-channel yielding an enhancement of the 
magnetic excitation intensity by a factor of 9 in the su- 
perconducting state as compared to the normal state |j2|. 
The difference between the in-plane and out-of-plane sus- 
ceptibility in the superconducting state arrises from the 
coherence- factor. 

Similar theoretical framework is currently used to de- 
scribe the spin excitation spectrum in spin-singlet HTSC 
cuprate superconductors f3(| . Theoretical works show 
that the opening of a d-wave order gap together with 
the exchange interaction lead to the appearence of a sim- 
ilar resonant feature below 2A at the antiferronagnetic 
wavevector. These theories successfully account for the 
observation by INS of the magnetic resonance peak in 
the superconducting state. 

Using unpolarized INS, one measures the superpo- 
sition of the out-of-plane and in-plane components of 
the susceptibility (see Eq. ||) and both components are 
equally weighted, when performing the measurements 
Qi=0. Thus, the predicted resonance feature should be 
observable, at the value of ^2A, if one obtains an exper- 
imental arrangement which allows to study the inelastic 
magnetic signal in this energy range. Due to the almost 
linear decrease of x"(qi, ^) towards low energies, see Fig. 
7, and due to the higher required resolution which implies 
less neutron flux, these experiments are extremely time 
demanding. We have analyzed the magnetic excitations 
in the superconducting phase on the cold triple axis spec- 
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trometer IN 14 at the ILL using a two crystal assembly, 
the results are shown in Fig. ^|. 

The right part of Fig. ^ shows a scan across the incom- 
mensurate peak at u>= 4 meV, i.e. the range already stud- 
ied with thermal neutrons. This signal can be determined 
with little beam time. Performing the same scan at 1.6 
and 0.8 meV requires considerably more time but still 
exhibits a well defined signal which seems not to experi- 
ence any change in the superconducting state at T=0.35 
K. The results of constant-Q-scans at Q=(0.7,0.3,0) are 
shown in the left part of Fig. |[ As there is no visible 
difference between the results obtained at 0.35 and 0.80 
K (both in the SC state), we have added the two scans 
below T c in the lower left part of Fig. |^. Importantly 
and despite efforts to get rather high statistics, there is 
no change visible in the energy spectra above and below 
the superconducting transition in the energy range of the 
superconducting gap. The spin susceptibility is not mod- 
ified appreciably across the superconducting transition; 
our data even do not show any opening of a gap. One can 
describe the energy dependence presented in Fig. 9 with 
the single relaxor using the same fitted parameters as the 
low temperature data in Fig. ||.a. However, the detained 
shape of the spin suceptibility (Fig. ^|) does not exactly 
match such a simple linear behavior but rather seems to 
indicate some anomaly near 2meV which requires further 
experimental work. 
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There is actually little known about the value of the su- 
perconducting gap in Sr 2 Ru04 . Laube et al. J57j have 
reported an Andreev reflection study where the open- 
ing of the gap is clearly visible in an astonishingly large 
energy range. The quantitative analysis of the spectra 
is quite involved; assuming a p-wave order parameter 
Laube et al. obtain 2A=2.2 meV which may be com- 
pared to the value expected within BCS-theory 2A= 3.55 
k B T c = 4.97 K = 0.43 meV. Our data shows that there 
is no change in the excitation spectrum for energies well 
below the reported value of 2 A, but it has not been pos- 
sible to investigate the lowest energies due to the strong 
elastic incoherent signal. With further increased resolu- 
tion {k f = 1.3 A" 1 ) we have scanned the energy range 
0.3-0.7 meV at 0.80 K again without evidence for a reso- 
nant feature. Furthermore, the comparison of two scans 
at constant Q=(1,0,0) did not yield any difference below 
and above T c meaning no ferromagnetic spin susceptibil- 
ity enhancement. 

The theory presented by Morr et al. Q should be 
considered as being quite reliable in the case of Sr 2 Ru04, 
since the RPA approach to the magnetic excitations is so 
successful in the normal state, and since Sr2Ru04 ex- 
hibits well defined Fermi-liquid properties at tempera- 
tures below 25 K. Therefore, the data in Fig. || gives 
strong evidence against a simple p-wave order parame- 
ter in Sr2RuC>4 with a maximum value of the gap of 
the order of the reported value J57|. However, in the 
meanwhile there are several indications that the order 
parameter is more complex. The recent specific heat data 
on the highest quality single crystals [^4| points to the 
existence of line nodes in the gap function which were 
then shown to be aligned parallel to the a,b-plane (hor- 
izontal line nodes). Such line nodes were explained by 
Zhitomirsky and Rice through a proximity effect between 
the active 7-band and the more passive one-dimensional 
bands [^7j . A modulation of the gap function along the 
c-direction will wipe out the resonance predicted for the 
non-modulated p-wave gap, since for the Q-position an- 
alyzed, (0.7 0.3 0), the electron hole excitation involves 
parts of the Fermi surface which are fully, partially or not 
gapped at all. In this sense the absence of any tempera- 
ture dependence in the magnetic excitation spectrum is 
consistent with the presently most accepted shape of the 
gap function. Further theoretical as well as experimental 
studies are required to clarify the possibility of a resonant 
feature at other positions in (Q,w)-space. 

IV. CONCLUSIONS 

Using assemblies of several crystals of Sr2Ru04 we 
have analyzed the magnetic excitations by INS. The in- 



commensurate signal arising from the nesting between 
the one-dimensional bands shows an asymmetry which is 
well explained by the full RPA analysis as a contribution 
mainly from the 7-band. 

The energy dependence of the incommensurate signal 
varies with temperature and exhibits a general softening 
of the spectrum upon cooling. This behavior indicates 
that Si'2Ru04 is approaching the corresponding SDW in- 
stability at low temperature even though this compound 
is not at a quantum critical composition. This inter- 
pretation is confirmed by the fact that the generalized 
susceptibility exhibits some w/T-scaling only above ~30 
K, i.e. in the temperature range where also the transport 
properties indicate non- Fermi- liquid behavior. 

The analysis of magnetic excitations besides the nest- 
ing ones shows only minor contributions. There is 
some evidence for additional magnetic scattering closer 
to the zone-center but still not peaking at the zone- 
center. This interpretation gets support from the fact 
that similar scattering is observed in nearly ferromag- 
netic Cai.5Sro.5Ru04 and from various RPA calculations 
which find excitations mainly related to the 7-band in 
this q-range. 

The magnetic excitations in Sr2Ru04 may be com- 
pared to the distinct types of magnetic order which have 
been induced by substitution. The dominant excitations 
reflect the SDW reported to occur in Si^Rui-^Ti^C^ at 
small Ti concentrations. The less strong excitations sit- 
uated more closely to the zone-center and most likely 
related to the 7-band become enhanced through Ca- 
substitution which drives the system towards ferromag- 
netism, but only for rather high Ca-concentration. 

Upon cooling through the superconducting transition 
we do not observe any change in the magnetic excita- 
tion spectra, which - combined with recent calculations 
p2| - indicates that the order-parameter in Sr2Ru04 does 
not possess simple p-wave symmetry. These experimen- 
tal findings are still in agreement with a p-wave order 
parameter modulated by horizontal line nodes. 
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